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The Chandra LETGS high resolution X-ray spectrum of the 
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Abstract. We present the Chandra LETGS X-ray spectrum of the nearby (~ 60 pc) neutron star RXJ1856.5- 
3754. Detailed spectral analysis of the combined X-ray and optical data rules out the nonmagnetic neutron star 
atmosphere models with hydrogen, helium, iron and solar compositions. We also conclude that strongly magnetized 
atmosphere models are unable to represent the data. The data can be explained with a two-component blackbody 
model. The harder component with temperature of kT^ ~ 63 eV and a radius i?2b — 2.2 km of the emitting 
region well fits the X-ray data and can be interpreted as radiation from a hot region on the star's surface. 
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1. The neutron star RXJ1856. 5-3754 



It was suggested by W alter et al. ( 1996| ) and supported by 
Neuhauser et al (|l997|) that RXJ1856. 5-3754 (or RXJ1856 
for short) is an isolated neutron star (NS). It has a strong 
and non-variable X-ray flux ~ 1.5 x 10~^^ erg s^^ cm~^ 
and a soft spectrum with a blackbody temperature kT^ ~ 
57 eV. Thousan ds of such NSs were expec ted in RASS 
(e. g., Colpi et al. |l993| ; Madau & Blaes |l994| ), moving fast 
through the interstellar medium. Only a few such objects 
have been found, three of them show pulsations on time- 
sca le of seconds (Haberl et al. 1997 , 1999 ; Hambaryan et 
al. ^001[ ). Neuhauser & Triimper ( |1999[ ) argued that the 
number of isolated NSs expected in RASS was overesti- 
mated mainly due to unrealistic velocity distributions. 

Walter & Matthews ( |l997| ) and Neuhauser et 
al. (1998) found an optical counterpart for RXJ1856 
with y ~ 26 mag. This and the large proper motion of 
~ 0.33 mas yr-i (Walter ^001^ Neuhauser ^OOlD are addi- 



tional arguments that it is indeed an isolated NS. Walter 
(2001) also detected parallactic motion, determined the 
distance to the source d = 6llg pc and suggested that 
RXJ1856 could be the remnant of the original primary of 
the C Oph system. This implies a NS age of ~ 10® yr. 

First spectral modelling of RXJ1856 based on its 
iJO^^T" data was presented by Pavlov et al. (1996; here- 
P96), who showed that the optical/UV flux pre- 



after 

dieted by NS atmosphere models depends drastically on 
the surface chemical composition. More recently Pons et 
al. (2002; hereafter P02) have analyzed combined optical 



and X-ray data and concluded that it may have either a 
Fe or Si-ash atmosphere. However, the inferred NS radius 
i? « 6 km and mass M « 0.9 Mq are not allowed for any 
plausible equation of state of the NS inner matter. 

RXJ1856 was observed with the Chandra Low Energy 
Transmission Grating Spectrometer (LETGS; Brinkman 
et al. 200C). First preliminary results were presented by 
Burwitz et al. ( |200lD . Here we describe the LETGS data 
(§2), results of spectral (§3) and timing (§4) analysis and 
discuss implications on the nature of RXJ1856 (§5). 



2. LETGS data extraction 

RXJ1856 was observed on March 10, 2000 with the stan- 
dard LETGS (LETG + HRC-S) configuration in a 56.1 ks 
exposure. The LETGS spectrum of RXJ1856 was ex- 
tracted from the level 1.5 event file. The only pulse-height 
filtering applied to the data was the removal of photons 
with a pulse-height amplitude equal to 255 as other pulse- 
height filters do not reduce the background level evenly 
for all energies. We used the extraction region recom- 
mended in the Chandra Proposers' Observatory GuideQ 
[POG]) for extracting the source spectrum (see Fig. |l|). For 
the background, large regions above and below the source 
extraction area are selected. Towards longer wavelengths 
A > 48.4 A the background regions become wider in order 
to maintain the constant ratio of areas between source and 
background spectral bins (here the ratio is equal to 6). The 
measured dispersed source count rate is 290±6ks~"'^ in 
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Wavelength [A] 

Fig. 1. Regions used for extracting the source and back- 
ground dispersed spectra. The vertical red lines on either 
side of the gaps between the HRC-S detector plates in- 
dicate the the areas that are masked out avoid getting 
artificial spectral features due to the dithering across the 
plate edges (see POG). 

the 0.15-0.82 keV range (where the source spectrum pre- 
vails over the background) . The extracted source spectrum 
(binned in 686 spectral bins) and the most up-to-date ef- 
fective area tables for the 1st order (status of October 31, 
2000)0 were used for spectral fits. 

3. Spectral analysis 

To search for spectral lines, we used the method applied 
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by Pavlov et al. ( |2001| ) to the dispersed LETGS data on 
the Vela pulsar. We first binned the extracted source-plus- 
background and background spectra in 0.02 A bins. Then 
we grouped 10-20 sequent bins and estimated the devia- 
tion of source counts in each bin of a given group from the 
mean value in the group. This method revealed maximum 
deviation at a 2.9 a level. 

Since there exist the ROSAT PSPC data (- 20, 200 
counts collected in a 6.3 ks exposure), we first checked 
whether the LETGS and PSPC data yield consistent re- 
sults in spectral fits. The blackbody fits showed that 
the allowed domains of the fitting parameters do not 
formally overlap with each other: kT^ = 57 ± 3 eV, 
= (3.7 ± 0.7) deo km (dgo = d/60 pc) and nH,20 = 
n]i/{\{)^° cm^^) = 1.47 ± 0.25 (the uncertainties are 
given at a 3cr level) with = 1.5 for the PSPC data, 
and A:T~ = 63 ± 3 eV, = (2.2 ± 0.3) deo km and 
nH.20 = 1-03 ± 0.20 [xl = 1.0) for the LETGS data. The 
difference between the blackbody models given by the two 
fits (see Fig. ||) is most likely attributed to uncertainties 
in the calibration of LETGS and PSPC. 

The next step in analysing the LETGS spectrum is 
to apply NS atmosphere models. The results of P96 and 
P02 showed that the light element, H and He, nonmag- 
netic atmosphere models can be firmly ruled out because 
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Fig. 2. The best blackbody fit to the LETGS spectrum 
(red curve) with the blackbody model given by the best 
PSPC fit (blue curve). 



the models (i) yield too small distances to RXJ1856, 
d ~ (5 — 6) (i?/10 km) pc and (ii) overpredict the op- 
tical flux by a factor of 20-30. P02 found that heavier 
element atmosphere models can be reconciled with the 
PSPC and optical data at kT = 57.4 eV, R 6.0(i6okm 
for pure Fe, and kT = 58.7 eV, R — 6.0^60 km for 
Si-ash compositions, assuming the gravitational param- 
eter g, = [1 - Rg/RV/^ = 0.77 (i?g = 2GM/^ is the 
Schwarzschild radius)Q Varying resulted in negligi- 
ble changes in the inferred model parameters unless the 
strong absorption line at (unredshifted) -Etruc — 0.79 keV 
(due to absorption by Fexviii ions) penetrates in lower 
(redshifted) energies E = grEnuc < 0.6 keV and signifi- 
cantly affects the quality of the PSPC spectral fits. The 
effects of gravity on the emergent spectra are too small 
to affect the results of the spectral fits. The fit of the 
PSPC spectrum with Fe atmosphere models of Pavlov 
& Zavhn ( |2000| ) yields d = 54 ± 10 pc, fcT = 45 ± 5 eV 
at i? = 10 km and M = 1.4Af0, that are well consis- 
tent with those obtained by P02. Although the LETGS 
spectrum revealed no statistically significant features, it 
is worthwhile trying the heavy element atmosphere mod- 
els on these data. To fit the LETGS data, we fixed the 
distance at d = 60 pc and obtained best fitting param- 
eters nH,20 = 1-72, kT = 40.3 eV, R = 12.0^60 km at 
Al — 1.4Mo (see red curves in the left panels of Fig. ||). 
The model fits well at (redshifted) £' < 0.5 keV where 
the atmosphere spectrum shows no strong spectral fea- 
tures. However, the situation is drastically different at the 
higher energies because of the strong absorption line at 
E — 0.81£'tiuc — 0.64 keV (see above) of a characteris- 
tic width ~ 70 eV and a weaker absorption complex at 
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^ For the atmosphere model fits we give parameters as mea- 
sured at the NS surface, T = and R = g,R°° 
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Fig. 3. The LETGS spectrum fitted with the Fe {left) and solar mixture {right) atmosphere models (see text for 
details). 



E ~ 0.55 - 0.58 keV. We note that the LETGS spectrum 
shows a strong instrumental feature at E = 0.61— 0.63 keV 
(POG). Despite of a formally small value of xt = 1-2 ob- 
tained in this fit, the deviation between the model and 
LETGS spectra at E > 0.5 keV are large enough to regard 
the fit as unacceptable. Increasing (i.e. reducing the 
redshift) does not improve the fit because the absorption 
lines remain within the observed energy range (see an ex- 
ample of the spectral fit at M = 0.2Mq, R — 9.8km and 
g-c = 0.97 in the upper left panel of Fig. ^ [blue curve]). 

The same strong Fe features are expected to be present 
in the Si-ash model spectra as Fe composes 68% of the ash 
(P02). Therefore, the situation with the Si-ash model fit to 
the LETGS spectrum should be similar to the Fe case. We 
also applied atmosphere models with a solar composition 
containing 2% heavy elements (Grevesse & Noels 1993| ). 
The model spectra show numerous prominent spectral fea- 
tures (mainly due to absorption by Fe, C and O ions) in 
the whole observed energy range (see right panels in Fig. ^, 
and also Rajagopal & Romani [1996|) which are inconsis- 
tent with the LETGS data. 



4. Timing analysis 

For the timing analysis we used 11,606 counts extracted 
from a 2"-radius circle centered at the zero-order image 



(we did not use the dispersed data as they are strongly 
contaminated by background). We had to restrict the 
search for a periodical signal at the (upper) frequency 
/ = 40 Hz because of the time-tag problem (see POG), 
unrecovered in this observation. As a result, the HRC tim- 
ing accuracy in this observation was about 4-5 ms (vs. 
16 /xs planned). We applied the standard Zl (Rayleigh) 
test (Buccheri et al. 1985| ) as well as the met hod b ased on 
the Bayesi an statistics (Greg ory h Loredo 1996| ; Zavlin 
et al. 2000| ; Hambaryan et al. 2001 ). The frequency range 
/ = 10^3 _ 40 Hz {P = 25 ms - 10^ s) revealed no sta- 
tistically significant signals (maximum significance is at a 
2.1(7 level). Assuming a sinusoidal signal, we put an upper 
limit /p < 8% on its pulsed fraction. 



5. Discussion and summary 

The most striking result of our Chandra LETGS data 
analysis (extended to the optical range [P96; P02]) is 
the fact that a simple blackbody model yields a much 
better fit to the RXJ1856 spectrum than either light or 
heavy element nonmagnetic atmosphere models. Note that 
the heavy element models are inconsistent with the high- 
resolution LETG spectrum. On the other hand, nonmag- 
netic H/He models do not work either because they pre- 
dict to high optical fluxes. A remedy may be provided 
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if RXJ1856 has a strong magnetic field as indicated by 
a bow-sliock nebula seen in optical observations (van 
Kerkwijk & Kulkarni 2000). Similar to the case with the 
nonmagnetic models, the magnetized H atmosphere fits 
result in a too small distance and a too high optical flux 
(P96). Magnetized Fe atmosphere spectra at kT 60 eV 
(Rajagopal et al. 1997) reveal a lot of strong absorption 
features in the 0.1-1.0 keV range which should have been 
resolved with the LETGS spectral resolution. Hence, we 
conclude that the available NS atmosphere models cannot 
represent both the X-ray and optical data on RXJ1856. 

Regardless of the chemical composition, radiation 
emitted from the surface of a magnetized NS should ex- 
hibit spectral lines of the electron (at Ebc — 11.6i3i2 keV) 
and proton (at Ebp — 6.3i3i2 eV) cyclotron resonances 
{Bi2 is the strength of magnetic field in units of 10^^ G). 
The lack of significant features in the LETGS spectrum 
(0.15-0.82 keV) may appear to exclude the magnetic fields 
B ~ (1.3 - 7.0) X 10^° G and B ~ (0.2 - 1.3) x lO^"* G. But 
we note that these restrictions on B are not very strin- 
gent because these lines may be rather faint if the surface 
layers are only weakly ionized. Besides, nonuniformity of 
the magnetic field over the surface should lead to a strong 
smearing of the lines (Zavlin et al. 1995). 

Remarkably, the X-ray spectrum of the other isolated 
NS, RX J0720.4-3125, obtained with XMM-Newton also 
shows no significant spectral features and is well fitted 
wit h a blackbody model of kT{^ ~ 86 eV (Paerels et 
al. |200lD . 

P02 argued that a two-component blackbody model 
can be reconciled with both the optical and X-ray data on 
RXJ1856. In this model the hard component of kT°° ~ 
55 eV is emitted from a ^ 20% fractional area on the 
NS surface and fits the X-ray data, whereas the soft com- 
ponent of kT°° ~ 20 eV represents radiation from the 
cool surface and matches the optical data. The required 
non-uniform distribution of the surface temperature may 
be due to the strong dependance of thermal conductiv- 
ity of the NS crust on the magnetic field a.t B > 10^^ G 
(Greenstein & Hartke |l98|) . At temperatures kT < 86 eV 
and magnetic fields B > 10^^ G hydrogen (if present on the 
surface) can be in form of polyatomic molecules and/or a 
condensed liquid (Lai & Salpeter 1997). Although, to our 



knowledge, no reliable calculations have been done, one 
may speculate that such a condensed matter surface emits 
radiation close to the blackbody spectrum at a tempera- 
t ure c lose to that of the surface, as suggested by Pavlov 
(200C)^. Then, the two-component blackbody interpreta- 
tion may be considered as a simplified model of the ther- 
mal radiation from such a magnetized and relatively cool 
NS. 

New important information for elucidating the na- 
ture of RXJ1856 is expected to come from the forthcom- 
ing 50 ks XMM-Newton and 450 ks Chandra observations. 
They will yield an order of magnitude more counts in the 
dispersed spectra than in the data presented above. This 
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will allow one to take more advantage of the high spec- 
tral resolution. Additionally, the new data would enable 
an unambiguous detection of X-ray pulsations from the 
source even if the pulsed fraction is as low as 3%. Once 
the period is found in the new data, it could be traced back 
to the earlier ROSAT a.nd Chandra data. This would give 
estimates on the period derivative and, as consequence, 
magnetic field and age of RXJ1856. The latter, compared 
with the estimate derived from the optical observations, 
would shed more light on the nature of this enigmatic 
source. 
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